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ABSTRACT 

There are several different ways to represent space- 
craft attitude and its time rate of change. For spinning 
or momentum-biased spacecraft, (me particular repre- 
sentation has been put forward as a superior para- 
meterization for numerical integration. Markley has 
demonstrated that these new variables have fewer 
rapidly varying elements for spinning spacecraft than 
other commonly used representations and provide ad- 
vantages when integrating die equations of motion. 
The current work demonstrates how a Kalman filter 
can be devised to estimate the attitude using these new 
variables. 

The seven Markley variables are subject to one con- 
straint condition, making the error covariance matrix 
singular. The filter design presented here explicitly 
accounts for this constraint by using a six-component 
error state in the filter update step. The reduced dim- 
ension error state is unconstrained and its covariance 
matrix is nonsingular. 

1. INTRODUCTION 

Attitude estimation for a torque-free spinning space- 
craft reduces to a comparatively simple problem if the 
angular momentum vector lies along the major prin- 
cipal axis. In this case, die spin axis remains constant 
in both the body frame and the inertial frame, and 
mean values for the spin direction, period, and phase 
can be determined using a batch method. 

It is more difficult to estimate a time-dependent atti- 
tude, such as occurs during maneuvers or when wire 
booms or other appendages are vibrating. Even in the 
absence of torques, the angular momentum will nutate 
in the body frame if it is not parallel to a principal 
axis. This effect can be important in case of a failure 
in the nutation damper. 

For these reasons, a new attitude determination utility 
for spinning spacecraft is being developed by the 
Flight Dynamics Analysis Branch at the National 
Aeronautics and Space Administration (NASA) God- 
dard Space Flight Center. This paper describes the 
Kalman filter design proposed for this utility. 

Attitude estimation methods for non-spinning, three- 
axis stabilized spacecraft often make use of the Euler 


symmetric parameters (commonly called the quater- 
nion) to represent the body attitude with respect to an 
inertial reference frame [1]. The quaternion is a four- 
component, globally nonsingular attitude representa- 
tion. The state vector typically has seven components 
consisting of the quaternion plus three elements that 
provide bias corrections to the measured rates. How- 
ever, spinning spacecraft usually do not carry rate- 
sensing gyros, so the rotation rate vector itself must be 
estimated rather than just the biases. Rates are needed 
for time propagation since the attitude dynamics 
satisfy a 2 nd -order differential equation. 

One disadvantage of using the quaternion for spinning 
spacecraft is that all four components generally are 
rapidly varying. It is shown in [2] that alternative 
descriptions can be devised that have fewer rapidly 
varying elements, and thus, are easier to integrate nu- 
merically. These elements will be referred to here as 
Markley variables. The main goal of this paper is to 
present a design for an extended Kalman filter to esti- 
mate the attitude using these new variables. This filter 
is expected to have superior convergence properties 
compared to one using the quaternion and rates to rep- 
resent attitude dynamics. The Markley variables used 
in this paper form a set of seven parameters, defined 
in Section 2. Section 3 presents the Kalman filter 
design. 

One major complication is that the Markley variables 
are subject to a constraint condition. This constraint 
makes the state error covariance singular. However, 
only six parameters are required to represent the atti- 
tude dynamics (e.g., roll, pitch, yaw, and their time- 
derivatives). It is always possible to define a six- 
element error state that is unconstrained. Section 4 
presents a new formulation of the filter based on a 
reduced, unconstrained, six-component representation 
for the state update and error covariance. Section 5 
gives a summary and plans for future work. 

2. MARKLEY VARIABLES 

The seven variables of the attitude dynamics repre- 
sentation described in [2] are the angular momentum 
in an inertial attitude reference frame (typically the 
geocentric inertial frame), I/, the angular momentum 
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in the body frame, L ft and a rotation angle, £ defined 
below. (Vectors will be denoted with bold characters.) 
These are subject to the constraint that the magnitude 
of the angular momentum vector is the same in the 
inertial and body frames. 


Let the A BE rotation define die angle £ so A BE can be 
written [1] 

4 * = cos £/ - ~^-[L b x] (7) 


The angle £ can be defined in a number of different 
ways. For this paper, Markley’s second representation 
will be used, as follows. An intermediate frame is 
defined by the rotation matrix 


A b m L 2 [L B ■ Ljl + thlflhlLL 

a L * ' l 2 + l b -l, 

-LjLl+L^]] 


( 1 ) 


where / is die 3x3 identity, L is die angular momen- 
tum vector, I is its magnitude, and the subscripts I, B, 
and E refer to the inertial, body, and intermediate 
frames. This definition is nonsingular as long as the 
angular momentum is nonzero and the body frame is 
not rotated exactly 180 degrees from the inertial 
frame. (In practice, this geometric limitation can be 
circumvented by introducing extra rotations to re- 
define die inertial reference frame as needed.) 


where the cross-product matrix is defined as 



for any vector v. 

2.1 Equations of Motion 

Expressions for the time dependence of the Markley 
variables are needed to propagate the attitude between 
sensor updates in the sequential filter. The angular 
momentum in the inertial frame satisfies the equation 
of motion 



( 9 ) 


The definition of Aei uses die angular momentum ex- 
pressed in both inertial and body frames. Caution is 
needed because of this mixing of frames; in the con- 
text of Eqn. 1, the Lb and L/ cannot be thought of as 
representations of an abstract vector. 

The matrix A a has the property that 

A EJ L i =L B (2) 

The attitude matrix, A B u is die transformation from the 
inertial to die body frame; thus, 

A B jL/ — L b (3) 


where N/ is the total external torque. The correspond- 
ing equation of motion expressed in the body frame is 

^- = N B -w m *L B (10) 

where 

(o tiI = J~'(L b -hg) (11) 

is the body rotation rate vector, J is the spacecraft 
moment of inertia tensor, and h B is die angular mo- 
mentum relative to the body frame of any internal 
moving components. The external torque expressed in 
the body frame is 


The Aei and A B i matrices both transform L/ to L B , but 
differ by a rotation about L B . This distinction is seen 
as follows. Define matrix A BE in terms of Aei and A B i, 

Abe s A B iA b (4) 

then multiply L; by A B i and use Eqn. 4, 

AseA^Lj — AgjLj (5) 

With Eqns. 2 and 3, this becomes 

AbE^B ~ I"B (d) 

Thus, A BE represents a rotation about L b . 


N b =A bi N, (12) 

Finally, [2] shows that £satisfies the equation 

d£ zfc + L / )- 0 >, / + L- 2 (L g xL / )-(jV g+ 7V / )] (13) 
dt 1? +L b - L, 

3. KALMAN FILTER 

There are two main parts that need to be specified for 
a Kalman filter design: the propagation step and the 
update step for both the state vector and the state 
covariance matrix. These steps are discussed in Sec- 
tions 3.1 and 3.2. Other details about Kalman filtering 
for spacecraft applications can be found in [3], for 
example. 



3.1 Time Propagation 


With this change, the expression for F is 


Section 2.1 presented the equations of motion for the 
Markley variables. These equations can be integrated 
numerically to propagate die full seven-component 
state forward in time between sensor observations. In 
the test version of the filter currently being developed, 
the time integration is implemented using a simple 
2 nd -order accurate Euler method (time-centered), with 
file total torque expressed as a sum of the command 
torques and environmental torques from gravity grad- 
ient and residual magnetic dipole perturbations. The 
constraint on the norm of L can be approximately 
maintained by using the average of the norms of Lb 
and L[ at each step, or L r can be forced to have norm 
exactly equal to that of Lb, as seen below. 

The state error covariance matrix, P, is defined as the 
expectation value of the outer product of the state 
errors, P - E[SX (S’ 7 ). This is a 7x7 matrix. However, 
since the norms of Lb and Lj are constrained to be 
equal, P is of rank 6. Anticipating this problem and 
the resolution that will be given below in Section 4, 
the full state vector now is defined to be 


X = 




(14) 


where Lj is the unit vector in the direction of the 
inertial frame angular momentum. This state vector 
still has seven-components, but the constraint now is 
that the first three elements have unit norm. 

The 7x7 error covariance, P, satisfies the equation 

— = FP + PF T + Q(t) (15) 

dt 

where (3(0 is the process noise, and F is obtained 
from the state dynamics equation, 

^ <i«> 

at 


with 



(tj ■ Njl + LjNj}/L 


^3x1 

F = 

®3x3 

\L b x]/ — \j L B x] 

^3x1 


d,L/(L 2 + Lj L b ) 

d B LHfi +Lj-L b ) 

^lxl 




(19) 

with 




d, 



(20) 

and 




d B = 

:-(L b+ L i J^ + {J-' 

l L B J-(N B + Nj[i B 

x]/L 


+(L B +Ljr l +(i r j- l i B )L T B 

-\f. B xL I \{N B +N l )£ B IL (21) 

3.2 Sensitivity Matrix 

The sensitivity, H, is the matrix of partials of the 
sensor observation with respect to the state X. Only 
vector observations are considered here. For example, 
tile output from V-slit Sun and star sensors and 
magnetometers are vectors. For these sensors, H is a 
3x7 matrix. 

Body frame vector observations can be modeled as 
v t = a bi v i + « 

= exp([«x]) 4 *v / +n (22) 

*(/+[ax])v fi 

= V B~[ V B *]« 

where v/ and v B are the reference vector expressed in 
the inertial frame and body frame, respectively, A at is 
the current estimate of the attitude matrix, a is the 
negative of the unknown attitude error vector, and n 
represents random noise with covariance R. Thus, the 
partial of v B obs with respect to a is -[vgx]. Using a in 
the error state rather than -a simplifies signs else- 
where when constructing a quaternion filter. 


dX 


Similarly, the negative attitude error, a, can be related 
(17) to the attitude quaternion, so that 


The function f(X,t) is given by Eqns. 9, 10, and 13, - = [g x] (23) 

except that Eqn. 9 now becomes da 2 2 



where q is the vector part of the quaternion and q 4 is 
(18) its fourth component. 



Next, the attitude quaternion can be expressed as a 
function of the Marldey variables, and the partials of q 
with respect to X can be determined. This yields 

dq C ; cos(4V2) f C 2 sin(C/2) (24) 

dL, L^2(L 2 +L r L B ) 

dq _ C cos(C/2) [ C 4 sin(C/2) (25) 

SL s L^2 {l 2 + L, ■ L b ) ^2{f}+L r L B ) 

and 


4.1 Unit Vector Covariance 

The 3x3 block of P corresponding to the uncertainty 
in the unit vector angular momentum can be related to 
die covariance corresponding to the 3-vector Lj as 
follows. If the estimated angular momentum and its 
error are written 

L’ = L + SL (32) 

where SL has zero mean, then the corresponding unit 
vector and its error are 


(£js x £,)sin( 472 ) . (/.„ + /., )cos(4~/2) ( 2 6) 
H 2L-yj2(L 2 +L r L b ) 2^2{^+L,-L b ) 


L' = (L + SL^L 1 +28L-L + SL 1 

* L+(l-LL r )^ 


(33) 


where 



where the subscript / has been dropped here 
simplicity. If the angular momentum covariance is 

for 


(27) 

p L = e[slsl t J 

(34) 

(L b +L,)LL t b / 2 

(28) 

then the covariance of its unit vector is 


‘ L‘ + L r L, 



p l = l 2 (i-ll t )p l (i-lL t ) 

(35) 

C (£ *1 {L’ 

r2 + t2 . r r 

l-M 

(29) 

4.2 Reduction to 6-Component State 



L 2 + L r L B J 


and 

r - 1 + (30) 

4 l 2 +l,-l b 

The 3x7 sensitivity matrix is obtained by combining 
all these pieces using the chain rule, 


Next, construct the matrix M that rotates the Z-axis of 
die inertial frame to be parallel to Lj. With this trans- 
formation, it is possible to reduce the dimensionality 
of the state simply by discarding the third component 
of M r L , , which is manifestly equal to unity. Only the 
X- and T-components are kept since they can deviate 
from zero in the 1* order at die update step. The state 
update is shown explicitly in Section 4.3. 


H = 


8< -[v xl! 

(dqX 1 

dq_ 

dq 

8q~ 

ax J 

J 

dLj ’ 

dL B ' 

’ d £. 


(31) 


To this end, define 

M s/cos5 + (l-cos^)n» r +sin5[nx] (36) 


4. REDUCED REPRESENTATION 

The state error covariance is a 7x7 matrix of rank 6. It 
can be difficult numerically to maintain P as a rank 6 
matrix during the filter propagation and update steps. 
To avoid this problem entirely, it is preferable to cast 
the filter in terms of a six-component error state, 
reducing the dimensionality of the covariance matrix 
and the filter update. 

Sections 4.1 mid 4.2 derive several relationships that 
are needed in Section 4.3 where the reduced form of 
the filter is presented. In particular, the reduced state 
vector is presented in Section 4.2. 


where 3 is the angle between the Z-axis and L h and 

(37) 


. Zxl. 
n - 


ZxLj 


With this definition, one has 

* A. 

MZ = L, 

by construction. Also, define 
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(38) 


( 39 ) 



One can show that 

T * *T 

Mm n m 2} M =I-L,L, 
Now, define die six-component state to be 


(40) 


JC = 


m n M T Lj 


sS T X 


(41) 


where S T is the 6*7 matrix 



m^U T 

®2x3 

®2xl 


^3x3 

I 

^3x1 


^1x3 

^1x3 

1 

_ 


(42) 


The m 2 3 matrix in Eqn. 41 explicitly removes the Z- 
component of the rotated unit angular momentum 
vector. Eqn. 41 also shows that 


dx _ &T 

w s 

One can show that S T S = / 6 , 6 and 
I M =SS r +Z L 


(43) 


(44) 


using Eqn. 40, and where / 6 „ 6 and / 7 * 7 are the 6^6 and 
7x7 identity matrices, and with Z L defined as 


is small over each short numerical integration interval. 
Recomputing S once for each time step should be suf- 
ficiently accurate for propagation of the covariance. 


Note that 


PZ L = Z L P = 0 


(48) 


since Lj is annihilated by the unit vector covariance, 
as shown by Eqn. 35. Thus, one can define 


f&s t fs 


(49) 


and 


Q*S t QS (50) 

to obtain a reduced covariance propagation equation. 


dP_ 

dt 


= FP+PF t + Q(t) 


(51) 


which is of the same standard form as Eqn. 15 and 
carries the same information content, but is not 
subject to any constraint condition. 


Similarly, the 7x3 Kalman gain 

k = ph t (hph t +r )~ 1 


(52) 


can be written in a reduced 6x3 form by introducing 
J 7 x 7 from Eqn. 44. Thus, 


lA 

^3x3 

®3xl 

^3x3 

®3x3 

^3x1 

^1x3 

^1x3 

®ixl 


4.3 Reduced Form for the Filter 


s t k=s t p(ss t +z l )h t 

[tf («SS r + Z l )p(sS t + Z l )h t + rY (53) 

= ph t [hph t +rY 


The 6x6 covariance matrix for the reduced state now 
can be defined as 


where R is die covariance of the sensor noise, n, in 
Eqn. 22, and where die reduced 3x6 sensitivity matrix 
is defined as 


P&S T PS 


(46) 


with time derivative 




dt dt 

~ s t f(ss t + z l )ps ( 47 > 

+ S T p(sS T + Z l )f t S + S T Q(t)S 


where die identity / 7x7 has been introduced in the form 
given in Eqn. 44. The time-dependence of S has been 
neglected. In die absence of external torques, S is 
exactly constant, and in practice, die time-dependence 


H&HS (54) 


Finally, die state update is 

Sx = K(v?-v?) (55) 

The 3 rd , 4 th , and 5* elements of the update vector Sx 
are added to the a priori estimate of L s . The 6 th 
element of Sx is added to the a priori estimate of £ 
The 1 st and 2 Bd elements of <5* are corrections to the 
unitized angular momentum vector in the rotated 
frame; the corresponding a priori elements of x are 



zero by construction. Thus, the updated inertial frame 
unitized angular momentum vector is 


L, =M 


5 : Cj 
Sx, 


( 1 m 


where M is given by Eqn. 36. The frill vector L t can 
be recovered at any time from L, and the norm of Lb- 


5. CONCLUSIONS 

An extended Kalman filter design has been presented 
for spinning spacecraft that makes use of the Markley 
variables. The Markley variables have much less vari- 
ability than other attitude representations such as the 
quaternion and rotation rate for spinning spacecraft, 
even in the presence of nutation or attitude man- 
euvers. This behavior makes diem more suitable for 
use in an attitude determination algorithm. 

The Markley variables are a seven-parameter set sub- 
ject to one constraint condition. The constraint leads 
to a singular error covariance matrix. A method was 
given for reducing the error state in the filter update 
step to six independent components, resulting in a 
nonsingular 6x6 error covariance matrix. 

The convergence properties of a Kalman filter for 
spinning spacecraft using die Markley variables are 
expected to be superior to those of a filter using the 
quaternion attitude representation. A prototype system 
is being developed to study die filter performance. 
Much additional work is needed before this filter can 
be migrated to the operational environment. In par- 
ticular, the filter must be able to accept the different 
measurement types used on spinning spacecraft. Only 
vector sensors were treated in Section 3.2. The form 
of the H matrix will be somewhat different for each 
sensor type; although, the steps leading to Eqn. 3 1 
will be similar for all sensors. 

Another feature that must be included in an oper- 
ational version of the filter is the ability to handle the 
180 degree restriction in die Markley variables. The 
denominator in Eqn. 1 goes to zero if Lb’Lj = -L 2 , that 
is, if the body frame is rotated 180 degrees from die 
inertial frame. It is possible to redefine the inertial 


reference frame with an extra rotation whenever the 
body-to-inertial frame angle approaches 180 degrees. 
Then, to reconstruct the usual attitude, it is only a 
matter of bookkeeping to keep track of the times these 
extra rotations were applied. A similar problem occurs 
in Eqn. 37 if Lj and Z are 0 or 180 degrees apart. The 
matrix M given by Eqn. 36 then is undefined. In these 
cases, M is either die identity or an arbitrary rotation 
by 180 degrees, and the component of L, to be dis- 
carded is already known to be parallel to Z. 

Another major enhancement would be to estimate 
sensor biases along with the attitude dynamics state. 
The ability to estimate these biases is already an 
important part of the current batch-method spinning 
spacecraft attitude determination system, where the 
spin vector is assumed to be constant for die entire 
batch. However, sensor biases may not be sufficiently 
observable in the time-dependent case. More work 
will be needed to study this possibility. Nonetheless, 
even die current filter design for estimating only the 
attitude dynamics will add a significant new capability 
to the attitude ground support system used for Flight 
Dynamics support for a number of spinning spacecraft 
missions at the NASA Goddard Space Flight Center. 
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